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Abstract: We study the optical force arising when isolated gold nanowire 
pairs and metamaterials with a gold nanowire pair in the unit cell are illu- 
minated with laser radiation. Firstly, we show that isolated nanowire pairs 
are subject to much stronger optical forces than nanospheres due to their 
stronger electric and magnetic dipole resonances. We also investigate the 
properties of the optical force as a function of the length of the nanowires 
and of the distance between the nanowires. Secondly, we study the optical 
force in a metamaterial that consists of a periodic array of nanowire pairs. 
We show that the ratio of the size of the unit cell to the length of the 
nanowires determines whether the electric dipole resonance leads to an 
attractive or a repulsive force, and we present the underlying physical 
mechanism for this effect. 
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1. Introduction 

The idea that electromagnetic radiation carries linear momentum goes back to James Clark 
Maxwell and his contemporaries Q); it is indeed the only way to reconcile Maxwell's equa- 
tions with the principle of conservation of linear momentum [2 |. Transfer of linear momentum 
from electromagnetic radiation to matter particles happens in the well-known phenomenon of 
radiation pressure, which is at least partially responsible for the tail of comets pointing away 
from the sun and is used in laser cooling to cool down atoms to temperatures close to the ab- 
solute zero [3-5 1. Recently, it has been proposed to harness optical forces, i.e., the forces that 
arise when linear momentum is transferred from photons to matter, in micro- and nanophotonic 
systems ^\U\- 



Optical forces are commonly classified as either gradient or scattering forces. The gradient 
force is an optical force that is perpendicular to the propagation direction of the excitation 
electromagnetic field and is well known for its use in optical tweezers, where strong laser 
beams generate a piconewton force that is used for the manipulation of small dielectric particles, 
including DNA, enzymes, and biological entities such as cells and bacteria. The physics of 
optical tweezers can be understood by recognizing that a dielectric particle can lower its energy 
by moving towards a region with higher field intensity [2|. More recently, researchers have 
exploited the gradient optical force for all-optical actuation of nanomechanical systems. The 
typical setup consists of two suspended waveguides close to each other; each of the waveguides 
sits in the exponential tail of the other waveguide's mode and will therefore be subjected to 
an optical force of a few piconewtons per milliwatt of optical power (6|[8j. The force can 
be enhanced using various types of optical resonators, such as ring resonators, up to a few 
nanonewtons per milliwatt 1 9 1 . 

As opposed to the gradient force, the scattering optical force imparts momentum parallel 
to the propagation direction of the excitation field. It has been studied extensively in the field 
of cavity optomechanics. The typical setup in this case consist of an optical cavity of which 
one mirror is free to oscillate as a mechanical harmonic oscillator. It has been demonstrated 
theoretically and experimentally that the optical force on the mirror and the resulting coupling 
between the optical modes and the mechanical motion leads to peculiar effects such as amplifi- 
cation and cooling of the mechanical oscillator 1 10 1. One of the goals of cavity optomechanics 
is to reduce the thermal occupancy of the mechanical oscillator such that the zero-point fluctu- 
ations of the — macroscopic — mechanical oscillator become observable. For our purposes, it is 
sufficient to note that the optical force can be strong enough to affect macroscopic masses. 

All the systems discussed above need objects that are larger than the wavelength of light 
in order to generate a measurable optical force. In a study exploring the use of the optical 
force for the aggregation of metallic particles into periodic arrays, Hallock, Redmond, and 
Brus have shown that optical forces can also be obtained when subwavelength metallic spheres 
are illuminated by laser radiation fTT) . When the laser frequency is chosen close to a collective 
plasmon resonance of the pair of metallic spheres, a much larger optical force can be excited 
due to plasmonic field enhancement 1 12 1. Early studies of optical forces on isolated nanowires 
include work by Rockstuhl and Herzig [ 13] and by Halterman et al. fl4) . 

In this article, we want to explore optical forces in metamaterial elements. In recent years, 
researchers have found that strong currents can be excited in arrays of subwavelength electro- 
magnetic resonators that are of quasistatic or plasmonic nature. Well-known examples of such 
metamaterial constituents are metal wires and split-ring resonators, i.e., metallic rings with a 
gap to provide sufficient capacitance. First fabricated in 2001 fl5) , metamaterials combining 
wires and split-ring resonators have been shown to produce media with negative permittivity 
and negative permittivity p6] - [T8| . Split-ring resonators and related structures are now well un- 
derstood and many studies have been devoted to their optimization 1 19-21 1. However, due to 
the saturation of the magnetic response above terahertz frequencies, they are not suitable for 
photonic systems. Here we will therefore use another prototype of a metamaterial element: 
metal nanowire pairs (see Fig.[T|. This structural element, which was independently proposed 
and fabricated by several groups [22-24], has been shown to exhibit a pronounced magnetic 
dipole resonance in the infrared and visible spectra. 

We will continue this article with a description of the methods used to calculate the optical 
force in the considered structures: finite-element simulations and the Maxwell stress tensor 
formalism for the calculation of the optical force. Thereafter, in Sec. [3] we will determine 
the optical force in an isolated nanowire pair illuminated by a laser beam. We investigate the 
strength of the force as a function of the geometric parameters of the nanowires and compare 
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Fig. 1. The nanowire pair that is considered in this article for the purpose of studying the 
optical force in metamaterials. The propagation direction and polarization of the incident 
electromagnetic wave is indicated. 



with the force between two spherical nanoparticles. Finally, in Sec. [4] we consider an array of 
nanowire pairs and we show that the interaction between neighbouring wire pairs can strongly 
modify the nature of the optical force, even if it influences the resonances of the individual wire 
pairs only slightly. 

2. Methods 

We use COMSOL Multiphysics to solve Maxwell's equations for the scattering problem of 
an incident plane wave illuminating the nanowire pair or the array of nanowires. COMSOL 
Multiphysics uses the finite element method to discretize Maxwell's equations. For the isolated 
nanowire pairs, we embed the structure in a large spherical simulation domain that is terminated 
with a perfectly matched layer. For the nanowire pair arrays, we use a cuboid simulation domain 
with periodic boundary conditions on the faces parallel to the propagation direction and ports 
at the faces perpendicular to the propagation direction. Gold is modelled as a linear medium 
with complex, frequency-dependent permittivity obtained from experimental data by Johnson 
and Christy |25) . 

After obtaining the full wave electromagnetic field solutions, we calculate the optical force 
on each of the nanowires using the Maxwell stress tensor formalism. The total force is the force 
on all charge carriers in the nanowire and is given by 

Fi = JJJ (pEi + e ljk JjB k ) dV, ( 1 ) 

v 

where e, ; ^ is the Levi-Civita symbol and V is a volume that fully encloses the nanowire. Since 
the metal is modelled by its complex permittivity, there are only bound charge carriers. Unfor- 
tunately, it is difficult (though not impossible) to implement Eq. ([T]l in COMSOL Multiphysics, 
because of surface charge discontinuities at the edges of the metal wires. It is, nevertheless, 
straightforward to transform Eq. ([T} into (2) 

*=-W&+JIT§| dV ' 

v v 



where S; = £qc EijkEjBk is Poynting's vector and Tjj is the Maxwell stress tensor defined by 



T {] = eo [EiEj - I SijE 2 ) + 1 (fi,^ - - 5 y B 2 J , (3) 

with 8ij the Kronecker delta. In the stationary harmonic regime employed in this article, the 
first term in the right-hand side of Eq. Q can be removed by averaging over an optical cycle; 
the second term can be transformed into a flux integral yielding 

(Fi) = ^(Tij)njdS. (4) 
s 

Extending the integration surface to a box that is larger than the nanowire, all fields are well- 
defined and the integration can be straightforwardly carried out with COMSOL Multiphysics. 

For the interpretation of the results, we will also make use of the electric and dipole moments 
of the current distributions in the nanowire pairs. We remind that all currents in the nanowires 
are modelled as polarization currents, so that both dipole moments can be obtained from the 
local electromagnetic fields available in COMSOL. The electric dipole moment is obtained 
from 

Pz = JJJ W = JJJ (D z - e Q E z ) dV. (5) 
v v 

The magnetic moment is determined from r x Jp, yielding 

m y = JJJ ( z ij7 ~ x §r ) dv = ~ ico JJJ [z(Dx ~ £()£x) ~ x (Dz ~ £o£z)] dy • (6) 

v \ / v 

3. Optical force in an isolated gold nanowire pair 

3.1. Physical nature of the optical force close to the resonances of a nanowire pair 

We start our study with an isolated gold nanowire pair with diameter of 25 nm, length of 100 nm 
and distance between the wires equal to 25 nm; such nanowires can be fabricated with aspect 
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Fig. 2. (a) The common and (b) the relative optical force exerted on a pair of nanowires 
by an incident field of 50mW/(im 2 . The length of the nanowires is L = 100 nm, their 
diameter is D = 25 nm, and the distance between the nanowires is d = 25 nm. The relative 
force exhibits two resonances: the magnetic dipole resonance with an attractive force and 
the electric dipole resonance with a repulsive force. 
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Fig. 3. Electric (a) and magnetic (b) dipole moment of the charge distribution in an isolated 
nanowire pair with length L = 100 nm, diameter D = 25 nm, and interwire distance D = 
25 nm. The insets show the charge and current distributions of the eigenmodes at 445 THz 
(electric dipole resonance) and at 410 THz (magnetic dipole resonance). 



ratios up to 7 by auto-catalytic growth from seeded metal salt solutions, identified by their par- 
ticle plasmon resonances |26], and self-assembled side-by-side using molecular linkers p7[ . 
The nanowire pair is illuminated with a plane wave with power density of 50 mW/ p.m 2 (see 
Fig.fllfor the direction of propagation and the polarization). Using the procedure described in 
Sec. [2] we calculate the forces exerted on both nanowires: F\ and Fz- For our purposes, it is 
more appropriate to consider the common force F com = (F\ +F2) acting on the center of mass 
of the nanowire pair and the relative force F K i — (F\ —F^jl from which we can determine 
whether the optical forces will increase the separation between the nanowires (repulsive force) 
or decrease the separation (attractive force). In Fig. |2|a)-(b), we plot both the common and 
relative forces for the 100 nm long nanowire pair. We see that the relative force displays two 
pronounced resonances: an attractive force around 410 THz and a repulsive force at 445 THz. 
We can understand the nature of these forces from Fig. [3ja)-(b), in which we plot the elec- 
tric dipole moment and the magnetic dipole moment of the charge distribution in the nanowire 
pair. At the resonance close to 410 THz, the magnetic dipole moment is maximal and we can 
therefore identify this resonance as the antisymmetric eigenmode. The corresponding charge 
distribution, sketched in the inset of Fig. [3jb), shows large opposite charges at the edges of the 
nanowires leading to an attractive Coulomb force. At the resonance close to 445 THz, the elec- 
tric dipole moment exhibits a maximum, and this resonance therefore corresponds to the sym- 
metric eigenmode with charges of equal sign at the edges of both nanowires [inset of Fig.[3|a)], 
resulting in a repulsive Coulomb force. Note that the magnetic resonance is much narrower than 
the electric resonance. This is due to the weaker radiation damping of the magnetic dipole. 

3.2. Optical force as a function of the nanowires' length and the interwire distance 

We continue our study with the question of how the strength of the optical force in nanowire 
pairs depends on the nanowires' length and the interwire distance. From the perspective of cir- 
cuit theory, we can already anticipate that the force will become stronger with longer nanowires. 
Indeed, the larger inductance of a longer nanowire pair yields a higher quality factor of the cir- 
cuit and, hence, larger charges at the ends of the nanowires. In addition, a longer nanowire pair 
will provide a stronger coupling to the incident field due to a stronger electromotive force. In 
order to test this hypothesis, we have plotted the average optical force exerted on nanowire pairs 
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Fig. 4. The relative optical force in an isolated nanowire pair as a function of length. All 
nanowires have diameter D = 25 nm and interwire distance d = 25 nm. Note that the force 
in the longer rods is exponentially larger than the force in the shortest sample (which has 
L — D = 25 nm and is therefore actually a pair of nanospheres). 



of different lengths in Fig. |4j we keep the nanowire diameter and the interwire distance fixed 
to 25 nm and the power density of the incident wave is still 50 mW/ [im 2 . We observe that the 
forces increase over several orders of magnitude when the length is increased from L = 25 nm 
to L = 125 nm. Note that a nanowire with L = 25 nm is actually a sphere, which allows to com- 
pare the optical force between metallic spheres to the force between nanowires. We conclude 
that the optical force is significantly stronger in nanowire pairs. In future work, it will there- 
fore be advantageous to make the nanowires as long as possible, of course within the limits of 
fabrication technology and the validity of the quasistatic theory employed in this article (i.e., 
the nanowires must remain smaller than the wavelength). Finally, we note that the attractive 
force disappears for nanorods shorter than 90 nm; this is because the electric and magnetic res- 




Fig. 5. (a) The relative optical force in an isolated nanowire pair as a function of fre- 
quency for different interwire distances. All nanowires have diameter D = 25 nm and length 
L = 100 nm. (b) The relative optical force at the resonance frequency as a function of the 
interwire distance. 



onances coincide for the chosen parameters and the electric dipole resonance dominates the 
charge distribution. 

For completeness, we have also studied the optical force as a function of the interwire dis- 
tance; the results are plotted in Figs. |5Ja)-(b). We observe that the overall peak value of the 
relative force increases with smaller interwire distance. This is in accordance with the explana- 
tion of the force in terms of the Coulomb force between the charges at the ends of the nanowires 
as explained above. The reader will note that the resonances frequencies shift, but this can again 
be readily explained with quasistatic circuit theory and is not the subject of this article. 

4. Optical force in a nanowire pair metamaterial 

We are now in good position to study a periodic array of nanowire pairs (see Fig. [6] for an 
illustration). The effect of stacking metamaterials constituents in periodic arrays has been stud- 
ied extensively in recent years. For most metamaterials, the coupling between neighbouring 
resonant element is merely renormalizing the circuit parameters (inductance, capacitance, and 
resistance) and although this coupling needs to be considered for the accurate design of meta- 
materials, the nature of the resonances remains unaltered [28 1. Only when resonant elements are 
strongly coupled by placing them extremely close or under special conditions such as coupling 
to dark elements, hybridized resonances with novel behaviour can emerge, eventually leading 
to effects such as electromagnetically induced transparency p9||34| . Such extreme coupling is 
not considered here and one could therefore expect that the optical force is not significantly 
modified by the stacking of the nanowire pairs in a periodic array. Surprisingly, our simulations 
show that this is not true. 

In Fig. [TJa)-(b), we plot the common and the relative optical forces on a 100 nm long 
nanowire pair with D = 25 nm and d = 25 nm. The results look similar to the forces obtained 
for the isolated nanowire pair. First, we check that the nature of the resonances has remained; 




Fig. 6. Schematic of the array of nanowires studied in Sec. [4] 




Frequency / THz Frequency / THz 



Fig. 7. (a) The common and (b) the relative optical force exerted on a pair of nanowires 
in a 2D periodic array with unit cell 100 nm x 200 nm when illuminated by a laser field 
of 50 mW/pm 2 . The length of the nanowires is L = 100 nm, their diameter is D = 25 nm, 
and the distance between the nanowires is d = 25 nm. The relative force exhibits two res- 
onances: the magnetic dipole resonance with an attractive force and the electric dipole 
resonance with a repulsive force. 



for the periodic array this can again be achieved by calculating the electric and magnetic dipole 
moments of the charge distribution in a single unit cell; the resulting dipole moments are shown 
in Fig.[HJa)-(b). Although the stacking of the nanowire pairs in a periodic array does have some 
influence on its frequency response, the main features remain: we observe the magnetic reso- 
nance at 410 THz, so that the attractive force (F re i < 0) at that frequency can still be attributed 
to the antisymmetric charge distribution associated with the magnetic dipole resonance as ex- 
plained above for the isolated nanowire pair. Similarly, the electric dipole moment exhibits a 
resonance at 450 THz, so that the repulsive force around that frequency can be attributed to the 
symmetric charge distribution associated with electric dipole resonance. 

We now determine the relative optical force on a nanowire pair in metamaterials with sue- 
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Fig. 8. Electric (a) and magnetic (b) dipole moment of the charge distribution in a 2D 
periodic nanowire pair metamaterial with unit cell 100 nm x 200 nm. The length of the 
nanowires is L = 100 nm, their diameter D = 25 nm, and the interwire distance equals 
D = 25 nm. 
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Fig. 9. The relative optical force in a periodic nanowire pair metamaterial as a function 
of the unit cell size height h (for smaller unit cell height, the ends of the nanowires in 
neighbouring cells come closer). All nanowires have length L = 100 nm, diameter D = 
25 nm, and interwire distance d = 25 nm. The unit cell size is 100 nm x h. Note that the 
optical force decreases significantly for metamaterials with smaller unit cell height. 



cessively smaller unit cell height (the dimensions of the nanowire pairs and the width of the 
unit cell remain unaltered), so that the ends of the nanowires in neighbouring unit cells come 
closer; the results are displayed in Fig. [9] We start by noting that the positions of the resonances 
shift to lower frequencies. This behavior can be explained with quasistatic circuit theory as, for 
example, in Ref. J28) : the effective capacitance associated with both resonances is increased 
by the capacitance between two collinear nanowires in neighbouring unit cells, resulting in a 
lower resonance frequency. 

More interestingly, however, is the strength of the optical force. We observe that the strength 
of the repulsive force is significantly diminished when the nanowire pairs are stacked in a 
periodic structure with smaller lattice constant. This cannot be explained based on our previous 
considerations, since circuit theory does not predict that the strength of the resonance would be 
changed so dramatically; this is confirmed by the electric and magnetic dipole moments (not 
shown here) that remain approximately constant at the resonance frequency. We can, however, 
understand what happens from the current distributions. If the distance between nanowires in 
adjacent unit cells becomes smaller than the distance between the nanowires of the same pair, 
the charges at the ends of two nanowires in adjacent unit cells start to cancel out and the charge 
distribution at the end of the nanowires will rather have the character of an electric dipole; the 
Coulomb force between the charges at nanowires' ends therefore start to satisfy an 1/r 3 law 
typical for electric dipoles rather than the inverse square law found before and the resulting 
Coulomb force is significantly reduced. 

For periodic arrays with larger nanowires, we find that the force associated with the mag- 
netic dipole resonance can even be reversed from attractive to repulsive. In Fig. 10 we display 
the relative force for a periodic array of nanowire pairs with length L — 200 nm, diameter 
D = 50 nm, and interwire distance d = 50 nm. It can be shown from the dipole moments (in- 
dicated by arrows in Fig. 10 1 that the large negative peak in the relative force for the samples 
with the largest unit cells (h — 480 nm and h = 360 nm) corresponds to the magnetic dipole 
resonance. For the larger unit cells, the force is therefore attractive in agreement with the phys- 
ical picture explained above. However, when we decrease the unit cell height, so that the ends 
of nanowires in neighbouring unit cells come closer, the force is first decreasing and finally 
reverses to become repulsive. It can again be shown from the dipole moments (indicated by 
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Fig. 10. The relative optical force in a periodic nanowire pair metamaterial as a function 
of the unit cell size height h (for smaller unit cell height, the ends of the nanowires in 
neighbouring cells come closer). All nanowires have length L = 200 nm, diameter D = 
50 nm, and interwire distance d = 50 nm. The unit cell size is 200 nm x h. Note that the 
optical force reverses from attractive to repulsive for metamaterials with smaller unit cell 
height. The arrows mark the position of the magnetic dipole resonance. 



arrows in Fig. 10 1 that this positive peak is still associated with the magnetic dipole resonance, 
and not with the electric dipole resonance as one would expect based on the repulsive nature of 
the force. We repeat that the Coulomb force between the charges at the ends of the nanowires 
in one pair turns into a weaker dipole-dipole force when the unit cell height is decreased. In 



order to explain the reversion of the force from attractive to repulsive (see Fig. 11 for an il- 
lustration)), we have to take into consideration the Lorentz force between the currents flowing 
in both nanowires. This force is repulsive, since the currents in both wires flow in opposite 
directions. The repulsive Lorentz force was negligible compared to the Coulomb force in all 
previously discussed samples, but now becomes first comparable and finally stronger than the 
attractive Coulomb force. 

It is important that the reader keeps in mind that the strength of the underlying resonance 
remains approximately constant, i.e., that the magnetic (or electric) dipole moment depends 
only slightly on the unit cell size. The examples above therefore show that the mere action of 
stacking a metamaterial element in a periodic array can dramatically change the optical force, 
even though the underlying quasistatic resonances are not altered significantly. 

5. Conclusions 

We have presented a study of the optical force in metamaterials. In the first part of this article, 
we have shown that appreciable optical forces can be generated in gold nanowire pairs when 
illuminated with laser radiation. Specifically, we have shown that gold nanowire pairs allow for 



optical forces that are significantly larger than those observed in pairs of nanospheres 1 1 1 12 1. 
We have presented a complete physical picture explaining the forces as Coulomb forces be- 
tween the charges accumulated at the ends of the nanowires. Repulsive forces originate from 
the electric dipole resonance of the nanowire pair and attractive forces originate from the mag- 
netic dipole resonance. 

In the second part of the article, we have considered a metamaterial consisting of a periodic 
array of nanowire pairs. As long as the unit cell size is large enough — and the nanowire pairs 
in neighbouring unit cells are well separated — we recover the same behavior as for the isolated 
nanowire pairs. However, with closer stacking of the nanowire pairs, we find that the nature of 
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Fig. 11. Charge and current distribution in the nanowire pair metamaterial. (a) If the height 
of the unit cell is substantially larger than the length of the nanowires, the force at the mag- 
netic dipole resonance frequency has the nature of a Coulomb force between the charges at 
the ends of the wires and is attractive for the magnetic dipole resonance, (b) If the height 
of the unit cell is only slightly larger than the length of the nanowires, the force at the mag- 
netic dipole resonance frequency has the nature of an electric dipole-dipole interaction; the 
repulsive Lorentz force between the currents flowing in opposite directions can now be of 
comparable size or even dominate the optical force. 



the optical force is altered dramatically, even though the underlying electric and magnetic dipole 
resonances are modified only slightly. In some nanowire pair geometries, it is even possible to 
turn the optical force associated with the magnetic dipole resonance from attractive to repulsive. 
Two effects contribute to such behavior: (i) As the distance between nanowires in neighbouring 
unit cells becomes smaller than the distance between nanowires in the same unit cell, the optical 
force changes from a Coulomb force between point charges into a much weaker dipole-dipole 
interaction, (ii) As the force between the charges at the ends of the nanowires becomes smaller, 
the Lorentz force between the currents in the nanowires will finally start to dominate, which 
explains the reversion from an attractive to a repulsive force. 
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